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We study the interaction of a massless quantized spinor field with the 
gravitational field of N parallel static cosmic strings by using a perturbative 
approach. We show that the presence of more than one cosmic string gives 
rise to an additional contribution to the energy density of vacuum fluctua- 
tions, thereby leading to a vacuum force of attraction between two parallel 
cosmic strings. 
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1 Introduction 



Cosmic strings predicted in the framework of various gauge theories with 
spontaneously broken symmetries could have been created at cosmological 
phase transitions in the early universe [p], |2|. Since cosmic strings are pro- 
duced at very large energy scales, one might expect a highly curved spacetime 
around them. However in the case of a static and straight-line cosmic string 
there exists a very simple exact solution of the Einstein field equations, that 
describes a locally flat conical spacetime around the string ||. Furthermore, 
it has been shown that one may also construct an appropriate exact solution 
of the Einstein equations for a snapping cosmic string, which serves as a 
source for spherical impulsive gravitational waves |4[]. 

Although the locally flat structure of the spacetime implies that straight- 
line cosmic strings will not exert any local gravitational force on surrounding 
particles, the particles "interact" with the global conical structure. It gives 
rise to the distinctive gravitational effects, such as lensing of distant objects, 
conical bremsstrahlung, etc. [||, ||. On the other hand, it is well known that 
the non-trivial topological structures restrict the modes of quantized fields 
propagating in locally flat spacetimes, thereby providing the appearance of 
vacuum boundary effects [[Jj] . One may, therefore, consider the conical space- 
time around static and straight-line cosmic strings as an attractive model for 
investigation the influence of gravitational field on the behaviour of quan- 
tized fields. In particular, such investigations for a massless quantized fields 
of different spins, propagating in the spacetime of a static cosmic string have 
been carried out in [8-10]. The propagation of a massless quantized scalar 
field in the spacetime of a snapping cosmic string was studied in [[ll[ . 

Recently in []12] we calculated the effects of vacuum fluctuations of a mass- 
less quantized electromagnetic field propagating in the spacetime of multiple 
cosmic strings represented by N parallel static (fixed), straight-line strings. It 
has been shown that the presence of more than one cosmic string provides an 
additional contribution to the energy density of vacuum fluctuations, which 
results in the Casimir-like force of attraction between two parallel cosmic 
strings. 

The aim of the present paper is to extend the calculations of [E2| to the 
case of a massless quantized spinor field. First, for an instructive purpose, 
we shall rederive the expression for the vacuum expectation value of the 
energy- momentum tensor of a massless spin- ~ field around a single cosmic 
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string. For the case of multiple cosmic strings we shall adopt a perturbative 
approach, in which the gravitational field of the strings is treated as small 
metric perturbations about the Minkowskian spacetime. We construct, at 
first-order in metric perturbations, the Feynman propagator for the spinor 
field and evaluate the energy density arising from vacuum fluctuations. We 
shall also show that at second-order perturbations of the metric, the effect 
of vacuum polarization of the spinor field gives rise to a force of attraction 
between two parallel static and straight-line cosmic strings. Throughout the 
paper we use geometrical units, in which G = c = 1 and % ~ 2.612 x 
10~ m cm 2 . 



2 A single cosmic string 

The metric of a static and straight-line cosmic string lying along the z-axis 
of the cylindrical coordinate system is given by the interval 

ds 2 = dt 2 - dz 2 - dr 2 - b 2 r 2 d6 2 (1) 

where 6=1 — 4/i, and \x is the linear mass density of the string. This 
metric is locally Minkowskian that is readily seen by passing to a new angular 
coordinate <p — > h9. We have 

ds 2 = dt 2 - dz 2 - dr 2 - r 2 dip 2 (2) 

However the global structure of this metric is a conical, as for the new az- 
imuthal angle ip we have the following range 

< (p < 2tt(1 - 4/i). 

The presence of such a "deficit" in the azimuthal angle provides boundary 
conditions for quantized fields in the metric (|2|), which result in Casimir-like 
distortions of the spectrum of vacuum fluctuations. The effect of vacuum 
fluctuations for a massless spin-| field has been calculated in || by using 
the method of Green's functions. Here we shall briefly reproduce this result 
using a different approach, involving the summation of the exact field modes. 

The propagation of a massless spin-| field in curved space-time is gov- 
erned by the Dirac equation 

i7 M V/ ^ = (3) 
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(4) 



where the 7 M matrices satisfy the anticommutation relations 

{7^,7"}= 2 <r. 

The spinor covariant derivative operator is defined as 

v„ = <9 M - r M 

where r M is the spinor connection. Let us consider a tetrad of vector field 
satisfying the relations 



Vab — 9\iv &a e b J 



e a e 6 



<r (5) 

For the sake of further convenience we shall use the Newman-Penrose tetrad 



(6) 



1 and m fl m* fl = — 1 . 



of null vectors [[13 

e£ = {P, n", m^, m*^} 

which satisfy the normalization conditions l^ l n fl - 
The symbol * denotes a complex conjugation. It is clear that with the choice 
of the tetrad (H the symmetric constant matrices rjab and r) ab have the only 
nonvanishing components 7712 = 7721 = —^34 = —Vu = 1, Vab = V ab - Since 
the ordinary constant Dirac matrices 7 a satisfy the condition 



{7 , 7 6 } = 2?T, 

then the coordinate-dependent matrices 7^ obeying the relation 
defined through the above introduced tetrad by 



r 



(7) 
will be 

(8) 



while the spinor connection Y ^ is defined as 



(9) 



where the semicolon denotes the covariant differantation with respect to the 
metric g^ u . The explicite form of the 7^ matrices written in terms of the null 
vectors is 

/ n» -m*"\ 
-m^ 
F m*^ 



V2 







(10) 



/ 
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Turning to the cosmic string metric (|2|) one may choose the following null 

vectors 



rrv 



' (0,0,-1,-' 



m 



1 
1 



;i, i, o, o) 



\/2 • r 
Let us now consider the four component spinor 



^(o,o, -1, - 



if, 



^2 
^3 



(12) 



then the Dirac equation (|3|) projected onto the null vectors (|TTJ) is decom- 
posed into the two independent pairs of equations 



and 



where 



and 



Dip! + (6*+ p)ifa = 

A^ 2 + (5+ P)iP! = 

A^ 3 - f3)fa= 

Di/j 4 - (5+ f3)ip 3 = 



(13) 



(14) 



* U V 



2V21 



D 
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Fdp = ^(dt~ d z ) A = n»dn = (d t + <9 2 
1 i 



(15) 



are the directional derivative operators. Combining the equations (|13f) and 
(|14D in an appropriate way we obtain the decoupled set of equations 

□s^s = (16) 
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where we have introduced the spin-weighted operator 
with the spin- weight s = ±| and the spin- weighted modes 

The equations fll6|) may be solved by the seperation of variables and the 
regular at r = solutions have the form 

e- J ^- fcz - m ^J|_ +s |(pr) (17) 

where J(pr) is a Bessel function and p 2 = u 2 — k 2 , v = 1/6. Since the 
transformation properties of spinors implies that a 2 7r rotation of a spinor 



changes its sign |14]| , the single- valuedness of the solutions (|17D requires that 
m = n — ~ , or m = n + | with n = 0, ±1 ± 2... . Having in hand these 
modes one can easily construct the complete sets of positive- and negative- 
frequency solutions u smkp and v smkp of the Dirac equation, which satisfy the 
normalization conditions 

fu s 'm'k' P '(x) 7 Usrakp^^gd 3, x = J v s > m > k > v > (x) 7° v s m kp (x) d 3 x 

= 6 SS , 6 mm , S(k - k') 5 (P~P') 

(18) 

where u = w 7 and v = v ' 7 are the adjoint spinors. 

The canonical quantization are performed in a usual way , by expand- 
ing the field operator i/j (x) on the complete sets of positive and negative 
frequency modes 

roc poo 

'<P( X )=J2 dk I dpp[a smk pU smkp (x) + bl mkp v smkp (x)] (19) 

s,m J -°° Jo 

where d smkp represents annihilation operator for particles and bl mkp is the 
creation operator for antiparticles. It should be noted that as the spacetime 
is flat everywhere outside cosmic strings one can define a vacuum state by 
choosing positive frequency modes with respect to the timelike Killing vector 
d/dt of the flat spacetime. 
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The energy-momentum tensor for spin- | has the form 

T„u(x) = ~ [?7(„V^ - (V(^) 7 ,)^] (20) 

The vacuum expectation value of this quantity can be evaluated by repre- 
senting it as a bilinear function of the fields and performing a renormalization 
procedure at the coincidence points x = x'. Substituting the expansion (|H^) 
into the equation (p0|) we obtain that its vacuum averaged < T 00 (x) > com- 
ponent are reduced to the form 

< T 00 (x) >= lim (d 2 - d t d t >)^r^ (21) 

x— *x ZOO 

where the two-points function U(x, x') is 

TT(t t'\ 1/ , roo r oo p -i\/k 2 +p 2 r+ik( 



Y j e lvm ^- ip) \ dpp dk 



2oo 8ir 2 7^ jo J-oo ^/k 2 + p 2 

[J\ m u+i\(pr)J\ mu+ i\{pr')+ J\ mu _±\{pr)J\ mv _^(pr')\ (22) 

and r = t — t' , ( = z — z'. The integrals over k and p in this expression are 
evaluated by means of the corresponding formulae given in [II]]. As a result 
we have 

U{X,X') V 1 _ -(^l)^,) / Wn _i )+ i| Wn _I)_|| 

(23) 

where 

r 2 + r' 2 + ( 2 -t 2 nr — 

u = — , 4 = u + Vu 2 - 1 

Z rr 

and we have taken m = n—\. The renormalization of the function U(x, x') 
is achieved by substracting from this expression its pure Minkowskian value, 
i.e at v — 1. Since the equation ([Zl]) involves only the derivatives over t and 
t' one can put (p — tp', z — z' and r = r'. Then the evaluation of the sum in 
is significantly simplified and we have 

g ^_ Wn _i )+ i| + rk(n _i } _i|^ = 2 tj ^Z±l!j (24) 

n=— oo ^ > 
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The further calculation of the renormalized quantity U(x, x') becomes straight- 
forward and substituting the result into the equation fl2lD we arrive at the 
following expression for the energy density of vacuum fluctuations 

< ^oo(x) >= (- 2 " 1) (7u 2 + IT) (25) 

This result is in agreement with that of given in |J. The other components 
of < T^v (x) > can be found by using the symmetry properties of the 
space-time ([TJ) along with the conditions = and T^K V = 0. 

3 Multiple cosmic strings 

If one makes a transformation of the radial coordinate r — > r\/b, then the 
metric ([I]) can be transformed into the the following form 

ds 2 = dt 2 - dz 2 - e~ 2A ^\dx 2 + dy 2 ) (26) 

where 

N 

A(x,y) = J2 / L^\nr l n = [{x - a,) 2 + (y - A) 2 ] 1/2 (27) 
i=i 

It turns out that this metric is the exact solution of the Einstein equations 



lq| , describing the space-time around N parallel static cosmic strings, pass- 
ing through the points Xi = (at$, Pi). The corresponding tetrad of null vectors 
for the metric (|26|) can be chosen as 

Z" = -j= (1, -1, 0, 0) n» = -±=(1,1,0,0) 

m» = -^=e A (0, 0, = ~ e A (0, 0, 1 - i) (28) 

Using these vectors in the equations ([8]) and (|9|) we find that 

7 °(x) = 7°, l\x) = 7 3 , 7 1 (x) = e A 7 \ 7 2 (x) = e A 7 2 , (29) 

where 7 °, 7 1 , 7 2 , 7 s are ordinary constant Dirac matrices, and the only 
non-vanishing components of the spinor connection are 

r\ = - i 7 V e A d y A r 2 = - l - 7 V e A d x A (30) 
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Now we shall evaluate the vacuum expectation value of the energy-momentum 
tensor fl20D in the metric fl2"6|). We start with the expression 



< T^x) >= -7 lim Tr bv(V„ -W v >)+ 7,(V M - V^)]S F (x, x') (31) 
where 

S F (x, x') = -i<0\T {^{x') ip(x)) |0 > (32) 
is the Feynman propagator which obeys the equation 

iT^fy - rjSpix, x') = 5\x - x') (33) 

The substitution of the relations (p9|) and (|30|) into this equation enables us 
to cast it in the form 

i-f a d a S F (x, x') = 5 A (x - x') + V S F {x, x') (34) 

where 7° once again denotes the flat space-time Dirac matrices, 

V = i[{\- e-^dA + (1 - e- A h a d a + \e~ V^A] (35) 

and the index A takes the values (0,3) = (t,z), while a = (1,2) = (x,y). 
In order to construct the solution of equation (|4|) we use a perturbative 
approach. For this purpose, we assume that the linear mass densities of 
the cosmic strings are sufficiently small, (fii << 1), which is indeed the 
case for realistic cosmic strings (fii ~ 10~ 6 ) Then the metric (^) can be 
expanded in powers of /ij, about a fixed flat background, and the potential 
(|35"D can be considered as a small perturbing term in the equation (B4l). This 



approximation allows us to write the solution of equation fl34|) as perturbation 
scries 

S F (x, x') = SP + SPVSP + SPVSPVSP +■■■ (36) 
The zeroth order free Feynman propagator is defined as 

s ™^W(0F e " ( " (37) 

where jk = ^ a k a . The successive terms in the expansion (^) correspond to 
the higher order in fii contributions to the free propagator. Expanding the 
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equation ( |35| ) in powers of /ij and taking it into account in we obtain 
that the first order corrections to the free Feynman propagator are given by 



fi£W) =6^5/ d*k5(k )5(k 3 )A(k)e-^ J -r^^ld 

[l A q A + 7ff " ^7*] 7(9 - fc) e-^" x,) (38) 

In order to evaluate the finite vacuum expectation value of the energy- 
momentum tensor (|3TD one needs to regularize it, for which we shall use the 
dimensional regularization procedure |lf|. The latter gives to the the van- 
ishing value for the zeroth order propagator 0, however, substituting the 
first-order propagator (|38"D into the expression (|3T]), for its < T^\x) > com- 
ponent at the coincidence limit x — > x' we find 



< T&\x) >= ^Tr / d*k5(k ) 5(h) A(k) e~ ik *' J 



q 2 (q — ky 



7% iq [ 1 A Qa + 7<? - Ipk ] j(q - k) (39) 

In the framework of the dimensional regularization procedure the integral 
over q is calculated, by performing an analytical continuation to d dimen- 
sional space and then the result is expanded about d = 4 — e, (e — > 0) [[14]]. 
Having done all these we obtain 

h N r 

< T${x) >= - ^^-l + 2 } ^> / ^ (A;2)1_f e ~ lkn (40) 

In obtaining of this expression we have used the explicit form for the Fourier 
component of the function A(x,y) 

\(k)= Stt^//,— - (41) 

i=i K 

It is important to stress that the expression fl4"D ) is finite, as the divergent 



terms are exactly cancelled by one another. Indeed, using the integral [12 



"*)V H ^"^^ (42) 
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in Eg. (pro|), in which d — 2 — e, and v = 1 — ~, we finally obtain the following 
result 

<T^)>=~T,^ t (43) 
i07r i=i u 

It is seen from this expression that at first order in fii, the contributions to 
the energy density of vacuum fluctuations of spin-| field around N parallel 
static cosmic strings are linearly summed. One can easily see that in the 
case of a single cosmic string the expression (^) coincides with fl2"o] ) taken 
at fi << 1. 



4 Vacuum force between two parallel cosmic 

strings 

Let us now proceed to the second-order metric contributions to the vac- 
uum expectation value of the energy-momentum tensor (|31|). In analogy 
with the case of electromagnetic fluctuations |12[] we shall show that the 



energy density of vacuum fluctuations of a massless spinor field involves a 
term which depends upon the seperation distance between cosmic strings, 
therefore, produces an attractive force between the strings. 

Using the expressions (|35|) andfl37|) in the expansion ( p6|) we find that 
the second-order metric corrections to the free Feynman propagator have the 
form 

sP(x,x) = f d*k5(k )5(k 3 )A(k)e-* kx ' [ d 4 P 5(p )5(p 3 ) A(p)e" 1 ^' 



256vr 8 

d*q 7? {l A< lA + iq - \lk) 7(g - k) 
q 2 (q — k) 2 (q — k — p) 2 

[l B q B +j{q-k-p)+ ]pp] !{q-k- p) e'^ x ~ x '> (44) 



Substituting this expression into the < T 00 > component of Eq. ( pi]) at the 
coincidence limit x = x', we calculate the traces of 7 matrices using the well- 
known theorems f\M, then taking into account the relation ([II]) we arrive at 
the expression 

•* TV N r j4l 
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/ f « S(p.) e-<-> / f {q _ kn f_ k _ p)1 "fa, *, P) 

(45) 

where 

A/"(g, fc, p) = 4(g A g A ) [4 q A q A + (g- £;) 2 + (k+p) 2 - fcp] + 2(g 2 - 2 fcg) 2g 2 ) 

- k 2 (2 g 2 + 2 fcg- p 2 - 2 A; 2 )- (Jfep) (g 2 - 3 k 2 + 2kq)+2 (pq) (q-k) (3q-k) 

We note that in this expression we are interested only in contributions to 
the vacuum energy density which are proportional to the products of the 
linear mass densities of different cosmic strings. It is clear that the latter 
describes the energy density of vacuum interaction between the strings. As 
for the remaining contributions, they form higher order corrections to the 
vacuum energy density given by fl43|). However first we need to regularize 
the expression fl45|) , for which we again use the dimensional regularization 
procedure. We evaluate the integral over q by performing an analytical con- 



tinuation to d = 4 — e dimensions using the scheme, described in [14|. After 
some straightforward algebra, we arrive at the following result 



fi n r d l k r d l n 

n2 



<^<*>>=-£E*E*/^-" ( ~'' / 



i=l N j=l 



f 1 dz x f 1 Zl cfe 2 [40r(-2+ J)5 2 ~f + dr(-l+ e -)B l -i 
Jo Jo 2 2 

+ C 2 r(|)S-f] (46) 
where we have introduced 

C x = -12Q(Q- k- p)+A(k 2 + p 2 + hp) 

o 

C 2 = Q 2 (-2Q 2 + 5kQ- k 2 -hp+ 3pQ) 

-kQ(k 2 + kp+ 2kQ + ApQ ) + \k 2 (p 2 + 2k 2 + 3kp + 2pQ ) 
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and 

Qa = k a Z 1 + (k a + p a ) Z 2 

B = k 2 z x {l -z x ) + {k + pfz 2 {\ -z 2 ) - 2k(k + p)zxz 2 

The energy of vacuum fluctuations per unit length of the strings may be 
evaluated by means of the formula 

E = J d 2 Xy/=g < T 00 (x) > (47) 

which at the second order metric perturbations takes the following form 

E = J dxdy[ < T ( 2) (x) > - 2A(x) < T$ (x) > ] (48) 

For the sake of certainty let us consider two parallel cosmic strings. Using the 
relations (flUD and ( fl6[) in the equation (|48|) we first carry out the integration 
over x and y, then the calculations of remaining integrals become straight- 
forward and keeping only the terms involving the product of the linear mass 
densities of the cosmic strings we find 

E mt = ^ vC~)j d 2 k {k 2 )~i e~ ika . (49) 

It should be stressed that this quantity is finite as the involved divergent 
at e — >• terms are compensated by one another. Indeed taking the double 
integration over k using the integral (|42D we obtain the expression 



Emt= (so) 

157T a z 

where a is the seperation distance between the cosmic strings. It is clear that 
the presence of this energy, will produce an attractive force per unit length 
of the strings given by 

F= _^^l (51) 
15tt a 3 V ; 

As we have already mentioned above the static and straight-line cosmic 

strings do not exert any local gravitational force on surrounding matter. Here 

we have shown that the propagation of a massless quantized spinor field in the 

spacetime of more than one cosmic string induces a force of attraction (|5T|) 
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between two cosmic strings, which falls of as the third power of the separation 
distance. The reason for this is the restriction of the modes of quantized field 
by the multiconical structure of the spacetime around the cosmic strings. 
Unlike the case of a single cosmic string, it is difficult to construct the exact 
modes of the field equations in the metric of multiple cosmic strings, so we 
have used a perturbative approach along with the dimensional regularization 



procedure. We note that the expression ([)!]) coincides with the corresponding 



result for a massless scalar field |Tj|, while 2 times smaller than the result 
for a massless vector field [12 |. 

It is important to stress that the above result is obtained within the one- 
loop approximation and therefore holds provided that the separation distance 
between the cosmic strings is much greater than the typical thicknesses of 
their cores. 
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